Descriptions of Carbon isotopes within relativistic Hartree-Fock-Bogoliubov theory 
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Within the relativistic Hartree-Fock-Bogoliubov (RHFB) theory, the structure properties of Carbon isotopes 
are systematically studied. In order to reproduce the experiment data, we take the finite-range Gogny DIS with 
a strength factor / as the pairing force. The self-consistent RHFB calculations with density-dependent meson- 
nucleon couplings indicate the single-neutron halo structures in both and "C, whereas the two-neutron halo 
in ^^C is not well supported. It is also found that close to the neutron drip line there exists distinct odd-even 
staggering on neutron radii, which is tightly related with the blocking effects and correspondingly the blocking 
effect plays a significant role in halo formation. 
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I. INTRODUCTION 

During the past decades, the radioactive ion beams (RIBs) 
have greatly extended our knowledge of nuclear physics, from 
which are obtained the critical data for nuclear physics, as- 
trophysics, as well as for testing the standard model. With 
worldwide and rapid development of RIB facilities, the in- 
vestigations of the nuclear systems under extreme conditions 
generate new frontiers in nuclear physics. For example, the 
exotic nuclei fV-'4\ have fascinated more and more interests 
due to the unexpected exotic modes therein. One of the rep- 
resentatives is the nuclear halo structure characterized by a 
dilute matter distribution contributed by several (in general 
two) loosely bound valence neutrons (or protons) surrounding 
a condensed core, which was first found in "Li I5j|. As the 
typical light nuclei, the Carbon isotopes have been devoted 
many efforts to probing the possible halo structure |@^8l and 
specifically recent measured reaction cross section of ^^C Q 
seems to assert a new two-neutron halo structure, which has 
also attracted fairly large interests from the community ifTOl - 

In fact, the exotic modes keeping found in the weakly 
bound nuclear systems also bring serious challenges on the 
reliability of the nuclear theoretical models. When extend- 
ing to the limit of stability of isotopes or isotones, the single 
neutron or proton separation energies become comparable to 
the pairing gap energy, such that the continuum effects can be 
easily involved by pairing correlations and play a significant 
role in determining the structure properties of exotic nuclei 
IfTSl lTB I. In terms of Bogoliubov quasi-particle, the relativis- 
tic Hartree-Bogoliubov (RHB) theory ifTsHTTll has unified the 
descriptions of relativistic Hartree (RH) mean field and pair- 
ing correlations, and consequently the continuum effects are 
involved automatically. Since the first self-consistent descrip- 
tion of nuclear halo structure in "LifTTl, the RHB theory has 
been successfully applied in predicting the giant halos in Ca 
|[T8l[l9| and Zr |[T4l[^|20l isotopes, as well as the restoration 
of relativistic symmetry flT] and superheavy magic structures 

ma. 

With the inclusion of Fock terms in the mean field, the 
relativistic Hartree-Fock-Bogoliubov (RHFB) theory with 



density-dependent meson-nucleon coupUngs Il23l provides a 
new self-consistent platform for the exploration of exotic nu- 
clei, e.g., predicting the giant halos in Cerium isotopes |19|. 
In addition, the inclusion of Fock terms has brought substan- 
tial improvements in the self-consistent description of nuclear 
shell structures |24 1 and the evolutions |,25. 26], the relativistic 
symmetry restorations Il24ll27ll28l . and the low-energy exci- 
tation modes Ii29l . 



In this work, the structure properties of Carbon isotopes, 
particularly the possible halo structures therein, will be stud- 
ied systematically within the RHFB and RHB theories. The 
contents are organized as follows. In the Sec. |llj we intro- 
duce the general formalism of the RHFB equations with finite 
range (Gogny) pairing force. In Sec. [Ill] the discussions are 
concentrated on the halo structures and odd-even staggering 
(OES) of Carbon isotopes. Finally, a brief summary and Per- 



spectives are given in Sec. IV 



n. 



THEORETICAL FRAMEWORK AND NUMERICAL 
DETAILS 



As generally recognized, the nucleon-nucleon interactions 
are mediated by the exchanges of mesons and photons. Based 
on that, the model Lagrangian contains the system degrees 
of freedom associated with the nucleon tfr, isoscalar scalar cr- 
meson, isoscalar vector w-meson, isovector vector p-meson, 
isovector pseudo-scalar ;7r-meson and the photon (A) fields 
[24. 30 1 . Following the standard variational procedure, one 
can get the equations of motion for nucleons, mesons, and 
photons, namely the Dirac, Klein-Gordon, and Proca equa- 
tions, as well as the continuity equation for energy -momentum 
tensor, from which is derived the system Hamiltonian. In the 
terms of the creation and annihilation operators (cl,,Ca) de- 
fined by the stationary solutions of the Dirac equation, the 
Hamiltonian operator can be expressed as 



^ = X clc/lTaf} c'lclc/^'Ca' ^ Vf^^^^,^„ 



(1) 



where Tap is the kinetic energy and the two-body terms 
y^^^,^,correspond with the meson- (or photon-) nucleon cou- 
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plings denoted by 0, 

Tap^ j dr{f,,(r)i-iyV + M)if,p(r), (2) 

Kfia'13' = J drdr'if/a(r)(fifj(r')r^{x,x') 

xD^{r,r')iff/s'(r')4ra'{r). (3) 

In above equations, r^(x, x') represent the interaction matrices 
associated with the cr-scalar, w-vector, p-vector, p-tensor, p- 
vector-tensor, 7r-pseudo-vector and photon-vector coupHngs, 
and D^{r, r') denotes relevant meson (photon) propagator, and 
M is the nucleon mass (for details see Refs. Il23ll24l[30ll ). 

Standing on the level of relativistic Hartree-Fock (RHF) ap- 
proach, the contributions from the negative energy states in the 
Hamiltonian (1) are neglected as usual, i.e., the so-called no- 
sea approximation |30|. The Hartree-Fock ground state |<l>o) is 
then determined and consequently is derived the energy func- 
tional E, i.e., the expectation of Hamiltonian with respect to 



l<l>o>=f]c!"|0). 



E =<(Do|//|Oo), 



(4) 



!=1 



where the index / denotes the positive energy states and |0) is 
the vacuum state. In the energy functional E, the two-body in- 
teractions lead to two types of contributions, i.e., the direct 
(Hartree) and exchange (Fock) terms. Within RHFB |23 1, the 
mean field part contains both types of the contributions, i.e., 
the RHF approach |Tll, whereas within RHB the Fock terms 
are neglected just for simplicity. 

For the open-shell nuclei, the pairing correlations, which 
lead to valence particles spreading over the orbits around the 
Fermi level, have to be taken into account. Different from 
simple BCS method |[32l . the Bogoliubov theory can unify the 
descriptions of mean field and pairing correlations in terms of 
Bogoliubov quasi-particles. It is of special significance in ex- 
ploring the nuclei far from the y6-stability line where the con- 
tinuum effects become essential and the simple BCS method 
may beak down. In the prior studies with both RHB and 
RHFB theories, it is already demonstrated that the scattering 
of the Cooper pairs into the continuum plays an essential role 
in the formation of the halo structures ifTSl [141 [T9l . 

Following the standard procedure of the Bogoliubov trans- 
formation l,33ii34J . the RHFB equation can be derived as. 



h{r,r')-A 
A(r, r') 



A(r, r') 
-h{r, r') + A 



= E 



Mr) 



(5) 



where ipu and ipv are the quasi-particle spinors and the chemi- 
cal potential A is introduced to keep the particle number on the 
average. For the single-particle Hamiltonian h{r,r'), it con- 
sists of three parts, i.e., the kinetic energy h^'", local potential 



h'^ and non-local one , 



(6a) 



h^ir, r') = [^Tir)y5 + ^or +P^sr]6(r - r'), (6b) 



h"{r, r') 



Yc(r,r')YF{r,r') 
Xc(r,r')XAr,r') 



(6c) 



Detials are referred to Refs. Il23l[30l . The pairing potential in 
the RHFB equation (|5]l reads as 



with the pairing tensor k 



(7) 



(8) 



For the pairing interaction V'', it is generally taken as a 
phenomenological form with great success in RHB theory 
||T6l[35| and conventional HFB theory |l36l[37l. In this work, 
we utilize the finite-range Gogny force DIS (381 with addi- 
tional strength factor / as the effective pairing interaction, 

V(r, r') =/ ^ e^^''-'">l'"'>\Wi + BiP'^ - HtP^ - MiPcrP^), 



1=1,2 



(9) 



where /i,, W,, B,, Hi, Mi(i = 1, 2) are the Gogny parameters. 

Due to the numerical difficulties originating from both RHF 
mean field and finite-range pairing interactions, the integro- 
differential RHFB equation (|5]l is solved by expanding the 
quasi-particle spinors on the Dirac Woods-Saxon (DWS) basis 
||39l , which can provide appropriate asymptotic behaviors for 
weakly bound nuclei. For the calculations of Carbon isotopes, 
the DWS basis parameters are taken as follows: the spherical 
box-size is fixed to 30 fm and consistently the numbers of ba- 
sis states with positive and negative energies are chosen as 48 
and 12, respectively. 



III. RESULTS AND DISCUSSION 

In this work, we calculated the Carbon isotopes from '"C to 
^^C by the RHFB theory with the effective interactions PKAl 
(241 PK02 |25|and PK03 |25|, and by the RHB theory with 
PKDD [40J and DD-ME2 L41J. For the odd Carbon isotopes, 
the blocking effects are taken into account. In practice, several 
orbits around the Fermi surface are blocked separately and 



TABLE I: Blocked neutron (v) orbits for odd Carbon isotopes in the 
calculations of PKAl, PK02 PK03, PKDD and DD-ME2. 
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FIG. 1: (Color online) Single- (5„: left panels) and two-neutron (52,,: right panels) separation energies for Carbon isotopes from '''C to -^C. 
The results are calculated by the density-dependent RHFB with PKAl, PK02 and PK03, and by density-dependent RHB with PKDD. The 
filled circles and open up-triangles denote the results calculated by taking the Gogny force DIS with/without prefix strength factor as the 
effective pairing interactions, respectively. As the references, the data extracted from Ref. 1431 are shown in filled squares. 



TABLE II: The mean square deviation from the data f43l for the binding, single- and two-neutron separation energies of the Carbon isotopes. 
The results are extracted from the calculations by PKAl, PK02 PK03, PKDD and DD-ME2 with different pairing strength factor /. 
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0.67 
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1.45 


0.68 


1.80 



the strongest binding leads to the ground state (42). In table 
[l] are shown the blocking configurations for the odd Carbon 
isotopes. 

Figure [TJpresents the single-neutron separation energies S„ 
(left panels) and two-neutron ones S 2,, (right panels) of Car- 
bon isotopes calculated with PKAl, PK02, PK03, and PKDD 
as compared to the data (in filled squares) [43 1. The results 
calculated by DD-ME2 are omitted because similar systemat- 
ics are found as PKDD. To get appropriate pairing effects, the 
comparison is performed between the results by taking the ef- 
fective pairing force as Gogny force DIS (open up-triangles) 
and the one with the additional strength factor / (filled cir- 
cles). From Fig. [Tjone can find that the modification on the 
pairing force brings some systematical improvements on both 
single- and two-neutron separation energies, especially in the 
vicinity of drip line. The optimized strength factors are deter- 
mined from the agreement with the data on the single-neutron 
separation energies as shown in Table [HI It is seen that the sys- 
tematics on single- and two-neutron separation energies are 
somewhat improved with the introduction of the strength fac- 
tor Specifically with the original Gogny DIS, PK02 can not 
reproduce the drip line ^^C which becomes bound with en- 
hanced pairing force (see Fig. [T]i. In Table [ill it ^Iso shown 
that the calculations with PK02 provide the best agreements 
with the data on both S,, and S2n, which may imply the most 



reliable systematics. Therefore we utilize PK02 as the rep- 
resentative to analyze the structure properties of Carbon iso- 
topes. 

Aiming at the possible halo structure in Carbon isotopes. 
Fig. [2|shows neutron and proton density distributions for even 
(a) and odd (b) Carbon isotopes. As shown in Fig. [2|a), it 
seems that the neutron densities tend to be more and more dif- 
fuse while not distinct enough to support the occurrence of 
halo structure, when close to the drip line. In fact the two- 
neutron separation energy of the ^^C is 1.56 MeV from the 
recent data [43,1, which suggests that the last two valence neu- 
trons are still bound too deep to spread over a fairly wide 
range. Hence ^^C may not be a good candidate of well- 
developed two-neutron halo structure. Whereas in Fig. l2](b) 
distinct evidence is presented for the halo occurrences iiP^C 
and ''^C, which represent much more diffused neutron distri- 
butions than ^^C. In fact, the strong evidence of the halo in '''C 
can be found from the parallel momentum distribution of '^C 
after the breakup of '''C [6[. As shown in Fig. [l]nearly zero 
single-neutron separation energies of '^C and '^C can be also 
treated as another evidence for the existence of single-neutron 
halo structure. While for ^'C the negative value of S„ leads 
to a diverged matter distribution, which might not be a bound 
nuclear system. 

To further illustrate the halo occurrence. Fig. [3ja) presents 
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FIG. 2: (Color online) Neutron and proton density distributions for 
even [Fig. (a)] and odd [Fig. (b)] Carbon isotopes. The results are 
calculated by density-dependent RHFB with PK02 and the pairing 
strength factor is adopted as / = 1.1. 



the contributions of the neutron density from different orbits. 
It is clearly shown that the dilute matter distribution at large 
radial distance is dominated by low- j state 2si/2 and the con- 
tinuum, in accordance with the evidence of halo occurrences 
in "Li flTl and Ca isotopes im. Consistently Fig. ^h) 
presents another direct evidence, i.e., the number of neutrons 
NR>r located beyond the sphere with radius r. From Fig. [3jb) 
it can be deemed that there exist evident single-neutron halo 
structures in '^C and '"^C due to fairly large amount of neu- 
trons spreading far beyond the neutron radii r„ . In contrast the 
values of A^s>,. in '^C and ^''C drop sharply with the increase 
of radius r, consistent with the neutron distributions shown 
in Fig. |2] Combining with the results in Fig. [3|a), one can 
find that both canonical state 2si/2 and the continuum present 
substantial contributions in the formation of halo, while dom- 
inated by the formal one due to its zero centrifugal barrier 

As the complemental evidence. Fig. |4] shows the neu- 
tron canonical single-particle energies for the Carbon isotopes 
from i^C to 22c. In Fig. gthe ultra thick bars denote the oc- 
cupation probabilities of each orbit. From Fig. |4] one can 
find that the valence orbits 2si/2 and liis/a are close to each 
another, which correspond to a fairly high level density and 
therefore induces strong pairing effects. Although both va- 
lence orbits are rather close to the continuum limit, the self- 
consistent RHFB and RHB calculations only support '^ '^C as 
the candidates of halo nucleus instead of even drip line isotope 
22c, which can be well interpreted by the blocking effects. 
In '^■'''C the neutron quasi-particle states 1^1/2 are blocked 
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FIG. 3: (Color online) (a) Contributions to neutron density (p„) 
from canonical neutron orbits (p"'^) and continuum states for "C, 
and (b) neutron numbers (A'k>, ) beyond the sphere with radius r for 
'^C£'**C£ "C and ^"C. The results are extracted from the calculations 
of density-dependent RHFB with PK02. r„ hq, r„ isq, r„ i9q and r„20c 
denote the root mean square neutron radii of '^C, '^^C'^C and ^"C, 
respectively. 



and the corresponding contributions are mainly mapped to 
the canonical single-particle orbit 2ii/2. Different from the 
even isotopes, the neutron staying on the canonical orbit 2si/2, 
mainly mapped from the odd-neutron on the quasi-particle 
state l5i/2, then becomes much less bound due to lacking 
the extra binding from pairing correlations which is also il- 
lustrated from S„ in Fig. [TJ^)- As a result, the probability 
density of the valence state 2s 1/2 tends to be much diffuser 
than those in even isotopes to develop the halo structure in 
'^•'"^C. In fact due to the blocking of s orbit, the continuum 
effects are also enhanced in odd isotopes because the neutrons 
occupied £^5/2 orbit can be only scattered into the continuum 
by pairing correlations. 

As we know, the pairing correlations play significant roles 
for the halo occurrences in the even nuclear systems, not only 
in stabilizing nuclide itself but also in developing the halos by 
scattering the cooper pairs into the low -lying sot p orbits. The 
typical examples are "Li, the drip fine isotopes of Ca, Zr and 
Ce. While for the even Carbons particularly 22c it seems that 
the extra binding from the pairing correlations makes s orbit 
too bound to get dilute matter distribution. On the contrary, 
the odd neutron in s orbit presents substantial contribution in 
the formation of halo structure, which can be seen from the 
odd-even staggering (OES) on the neutron radius. Figure |5] 
(a) shows the neutron root mean square radii from '"^C to 22c 
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FIG. 4: (Color online) Canonical neutron single-particle energies 
for Carbon isotopes. The results are extracted from the density- 
dependent RHFB calculations with PK02. The length of thick bars 
denotes the occupation probabilities of neutron orbits and Ef in short 
dashed lines represent the Fermi energies. 



extracted from the calculations with PKAl, PK02, PK03 and 
PKDD. It is clearly shown that all the selected effective in- 
teractions except PKAl present distinct OES on the neutron 
radius. Specifically, as seen from Fig. [5] (b), such OES is 
determined by the valence orbits 2s 1/2 and also depends on 
whether the corresponding Bogoliubov quasi-particle s orbit 
is blocked or not. In practice the orbits around the Fermi sur- 
face of the neighboring even isotopes are blocked separately 
and the blocking configuration is then determined as the one 
leading to the strongest binding, i.e., the ground state. From 
table |l] one can find that the blocking configurations are con- 
sistent with the OES in Fig. [s] In '^C which has larger neutron 
radius than '^C in the calculations with PKAl and PKDD, the 
quasi-particle s orbit is blocked. Due to similar reason and 
consistently with the halo occurrence, the neutron radii of halo 
nuclei '^ '^C are distinctly larger than the even neighbors. Dif- 
ferent from PK02, PK03 and PKDD, PKAl in fact does not 
suppoit the halo occuiTence for '^C in which (1^/2 is blocked 
for the odd neutron and due to the pairing effects the neutrons 
occupied in the low- j state, i.e., the s orbit, are bound too deep 
to distribute extensively. 

It is well known that pairing correlation plays an impor- 
tant role in stabilizing the finite nuclei, especially the exotic 
ones. For "Li, Ca, Zr, and Ce isotopes, it is already demon- 
strated that the pairing correlations show positive effects in 
both stabilizing and developing the halo structures. While in 
the RHFB (except PKAl) and RHB calculations of '^ '^C, the 
quasi-particle s orbit is blocked and the corresponding contri- 
butions are mainly mapped to the canonical orbit 2si/2, which 
plays the dominate role in the formation of single-neutron halo 
structure of '^ '''C. This implies that the unpaired odd neutron 
may also contribute to develop the halo structure when it is 
not so deeply bound. In fact from the analysis above it is just 
due to the lack of extra binding from the pairing correlations 
that the odd-neutron can spread over far beyond the center of 
nucleus. 




FIG. 5: (Color online) (a) Neutron root mean square radii calculated 
by density-dependent RHFB with PKAl, PK02 and PK03, and by 
density-dependent RHB with PKDD, and (b) contributions from the 
core orbits (Lsi/2, lp3/2 and lpi/2), valence orbits (2^1/2 and W5/2) 
and the continuum. The results are provided by the calculations of 
density-dependent RHFB with PK02. 



IV. SUMMARY AND PERSPECTIVES 



In this work we have systematically calculated the ground- 
state properties of Carbon isotopes using the density- 
dependent relativistic Hartree-Fock-Bogliubov (RHFB) the- 
ory with PKAl, PK02 and PK03 as well as the density- 
dependent relativistic Hartree-Bogliubov (RHB) theory with 
PKDD and DD-ME2. It is shown that neutrons in the 2si/2 
state in '^C and '^C are weakly bound and direct evidences 
are obtained for the single-neutron halo occurrences in the odd 
Carbon isotopes '^C and '^C. In addition, the self-consistent 
density-dependent RHFB calculations do not support two- 
neutron halo structure in ^^C, disagreeing with the indications 
of the experimental reaction cross section measurement |9|. 
Pairing correlation is one of the important stability mecha- 
nisms of finite nuclear systems, and the blocking effect ac- 
companying with it also plays a significant role in determining 
micro structure for odd-A nuclei. In this work, because of the 
blocking effect, the neutrons in the 2si/2 state in '^C and '^C 
are unpaired and weakly bound. Compared to the calculation 
with PKAl, one can see clearly that when the qusi-particle 
state lds/2 is blocked in '^C, there dose not exist halo struc- 
ture. In a word, blocking effect plays a crucial role in the halo 
occurence for the odd-A nuclei. 
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